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ABSTRACT 


An  emphasis  is  placed  on  the  real-time  correction  of  ref raction- induced 
tracking  errors.  With  radar  data  from  three  missile  trajectories  the  cor¬ 
rected  range  and  angle  data,  using  real-time  equations,  is  compared  with 
ray-traced  data.  It  is  shown  that  these  eouations  can  be  effectively  used  to 
reduce  the  range  and  angle  errors  to  about  the  resolution  of  the  tracking 
radars;  that  is,  approximately  two  feet  in  range  and  five  thousandths  of  a  # 

degree  in  angle.  It  is  also  shown  that  the  effect  of  an  elevated  layer  on  track¬ 
ing  angle  errors  can  be  compensated  using  these  equations. 
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SECTION  I 


INTRODUCTION 

In  a  previous  report1  an  analysis  was  made  of  the  refraction-induced 
errors  produced  by  time  and  space  variations  of  the  medium.  These  ap¬ 
parent  variations  limit  onefe  ability  to  define  a  single  refractivity  profile 
which  can  be  used  to  compensate  these  errors  in  near  -  real- time .  Within 
these  restrictions  a  radiosonde  profile  obtained  close  to  the  time  of  launch 
and  at  a  location  in  front  of  the  tracking  radar  provides  a  first-order  des¬ 
cription  of  the  propagation  conditions. 

To  determine  the  significance  of  the  refractivity  fine  structure  ob¬ 
tained  from  the  radiosonde  requires  ray-tracing  analysis.  Due  to  the  time 
required  to  obtain  ray-traced  data  for  real-time  corrections  one  must  store 
the  error  data  obtained  from  the  tracings  in  advance  and  then  reacquire  this 
data  or  generate  simple  empirically-determined  functions  which  can  yield 
the  equivalent  error  data  in  real-time.  The  ray-traced,  storage  approach 
requires  making  use  of  a  large  part  of  the  tracking  computer  memory  bank. 

If  several  radars  are  required  to  provide  raw  data  into  this  central  computer 
this  approach  generally  becomes  unacceptable. 

The  second  method  using  functions  fitted  to  the  ray-traced  error  data 
has  several  variations.  If  the  empirical  functions  are  derived  directly  from 
error  data  the  resulting  equations  can  lose  identification  with  the  real-world 
parameters  which  produce  refraction  errors,  such  as  the  effect  of  the  range, 
elevation  angle,  surface  value  of  refractivity,  etc  These  equations  are  then 
mathematical  expressions,  probably  polynomials,  but  can  be  made  sufficiently 
simple  to  be  used  for  real-time  corrections.  By  preparing  a  library  of  these 
expressions  for  a  wide  range  of  measured  refractivity  profiles  experience 
should  permit  one  to  select  functions  to  meet  most  situations. 

Another  approach  is  to  assume  that  on  the  average  the  refractivity 
profile  can  be  expressed  by  a  simple  exponential  function2.  With  suitable 
approximations  equations  can  be  developed  for  real-time  range  and  angle 
corrections.  To  date,  the  development  of  these  equations  does  not  permit 
one  to  include  an  elevated  layer  in  the  model. 

The  following  report  describes  the  development  and  application  of 
some  new  correction  equations  which  can  be  used  in  real-time.  These 
equations  represent  continuous  functions  which  maintain  a  physical  identifica¬ 
tion  with  the  parameters  which  effect  refraction- indue ed  errors.  Inherent  in 
these  equations  is  the  ability  to  include  the  effect  of  an  elevated  layer  thereby 
making  these  equations  particularly  useful  for  low-angle  tracking.  A  com¬ 
parison  is  presented  of  the  differences  between  tracking  parameters 
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determined  from  these  real-time  equations  and  the  AFWTR  post-flight  ray¬ 
tracing  analysis  for  two  vehicle  launches  and  with  three  AFWTR  radars 
tracking  each  launch.  In  addition,  one  comparison  is  made  with  the  GERTS 
real-time  correction  program.  However,  no  inversion  layer  was  present 
for  this  launch 

The  comparative  results  indicate  that  meaningful  corrections  to  re- 
fraction- indue ed  tracking  errors  can  be  obtained  with  the  SURC  real-time 
equations.  If  accurate  tracking  data  is  required  only  above  five  degrees 
elevation  angle  these  equations  require  only  a  single  value  of  surface  re- 
fractivity  which  can  be  obtained  from  the  locally  associated  radiosonde 
launch  station  To  improve  the  accuracy  below  five  degrees  elevation  angle 
the  characteristics  of  the  inversion  layer,  if  present,  need  to  be  included. 

These  c onsiderations  indicate  that  if  further  investigations  are 
necessary  they  should  be  directed  toward  the  measurement  of  the  residual 
errors  within  an  experiment  where  corrected  radar  data  is  confronted  with 
known  target  position  data.  Suitable  test  targets  could  be  established  by  first 
using  high  altitude  aircraft  equipped  with  accurate  radar  altimeters  for 
short  range  studies.  Finally,  using  satellite  targets  of  known  trajectory,  the 
long  range  errors  at  both  low  and  high  angles  could  be  studied.  However,  in 
the  light  of  presently  available  tracking  accuracies  using  the  FPS-16  radars 
it  is  believed  that  the  bulk  of  residual  errors  will  be  largely  due  to  the 
electromechanical  characteristics  of  the  radars  and  that  the  residual  re¬ 
fraction  errors  will  be  almost  impossible  to  detect  in  a  systematic  way. 
Therefore,  further  intensive  investigations  of  the  nature  of  all  residual 
errors  is  deemed  to  be  meaningless  until  the  tracking  system  can  resolve 
range  to  less  than  two  feet  and  angles  to  less  than  0.  1  milliradians . 
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SECTION  II 


ANALYSIS 


It  was  shown  that  for  propagation  in  a  simple  exponential  atmosphere 
the  range  error,  AR,  and  true  elevation  angle,  B0,  can  be  expressed  by 
(figure  1). 


AR  - 


10'6  N s 
c  sin  y 


[i 


g  .  -(k2  +  2kg)i 

k  +  g 


(i) 


R 


B0  —  sin 


-1 


[—  f  sin  (90  -  r)  dR] 


(2) 


where  the  magnitude  of  the  ray -bending,  T,  is 

T~  »°'“NS  .  _  _g_  -(k2  +  2kg), 

tan  y  k  +  g 

and 

,  ,  ^  r «  «  "ko R-1  —  k»  0  p, 

sin  y  =  sin  0O  +  [k0  +  k2  e  3  ]  e  10 


(3) 

(4) 


g  = 


tan  y 


(5) 


k  = 


cos  y 


(6) 


N  =  N  s 


e-c(h  -  Hs) 


(7) 


The  empirical  constants,  k0,  .  .  .  ,  k3,  were  required  at  short  ranges 
and  small  elevation  angles  to  provide  ray-bending  calculations  from  (3)  in 
agreement  with  the  CRPL  exponential  model  atmosphere.  Their  variation 
with  station  value  of  r ef rac tivity,  Ns,  is  shown  in  figure  2. 


Since  the  CRPL  model  was  developed  to  correct  ray-bending,  the 
magnitude  of  the  exponential  decay  constant,  c,  (7),  tends  to  increase  with 
Ns.  2  The  range  error,  £R,  being  the  integral  of  N  tends  to  be  directly  de¬ 
pendent  on  the  surface  value,  Ns.  When  the  station  value,  Ns,  exceeds 
about  340  units  the  effect  of  c,  (1),  causes  the  range  error,  £R,  to  decrease 
as  Ng  increases,  which  is  abnormal.  However,  the  error  produced  by  this 
effect  is  very  small  compared  with  the  position  error  produced  by  the  effect 
of  ray-bending.  Therefore,  it  is  meaningful  to  concentrate  on  the  reduction 
of  target  position  error  produced  by  ray-bending  and  for  this  reason  the 
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FIGURE  1 .  RAY  GEOMETRY 
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FIGURE  2.  VARIATION  OF  EMPIRICAL  CONSTANTS  WITH 
SURFACE  REFR ACTIVITY 


CRPL  model  was  adopted  in  this  development.  Furthermore,  the  velocity 
error  is  more  seriously  affected  by  ray-bending  than  by  the  range  error.3 

The  true  elevation  angle,  B0,  can  be  calculated  directly  with  a  digital 
computer  by  writing 

R  N 

f  sin  (0O  -  T)  -  F  sin  (0Q  -  T* )  •  6Rt  (8) 

o  i=  1 


where,  T^,  is  evaluated  at  the  mid-point  of  each  interval,  6R{,  using  (3)  and 


6R{ =  R/N 


(9) 


Since  the  bending,  T,  is  a  logarithmic  function  of  range,  R,  it  is 
possible  to  greatly  reduce  the  number  of  elements  in  the  summation  by 
expressing 


Ml-i- 


(10) 


where 

ARt  =  R  [1  -  logN  (N  -  i  +  1)] 


(ID 


1=2,  3, 


N 


For  example,  whereas  25  intervals  is  required  with  (8)  to  maintain  a 
comparative  angle  accuracy  of  0.05  milliradians ,  the  number  of  intervals 
can  be  reduced  to  9  (N  equals  10)  with  (10),  This  represents  a  considerable 
reduction  in  execution  time. 


The  above  equations  can  be  used  in  real-time  to  determine  tracking 
errors  from  measured  radar  parameters  Appendix  I  shows  how  the  effect 
on  ray-bending  produced  by  an  elevated  layer  can  be  included.  This  first- 
order  layer  correction  still  permits  these  SU  AC  equations  to  be  used  in 
real-time . 
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SECTION  III 


COMPARISON  OF  TRACKING  ERRORS  BETWEEN  SURC 

CALCULATIONS  AND  AFWTR  RAY-TRACED  RESULTS  - 

OPERATION  5265 

Figures  3  and  4  show  the  true  elevation  angle  and  range  from  the 
South  Vandenberg  tracker,  (FPS-16,  023002)  for  a  vehicle  launched  from  the 
Vandenberg  firing  site.  This  angle  and  range  data  was  obtained  by  correct¬ 
ing  the  raw  radar  data  using  the  AFWTR  (post-flight)  ray-tracing  program. 

It  should  be  noted  that  this  program  was  not  used  in  real-time.  The  propaga¬ 
tion  conditions  were  defined  by  a  radiosonde  launch  from  South  Vandenberg 
close  to  vehicle  launch  time.  Figure  5  shows  the  lower  part  of  the  refrac- 
tivity  profile,  and  it  is  apparent  that  no  effective  layer  was  present  above  the 
height  of  the  Tranquillan  Peak  radar.  The  radiosonde  data  for  all  associated 
radars  and  launches  is  included  in  Appendix  II. 

The  raw  radar  range,  R,  and  apparent  elevation  angle,  0O,  were  used 
as  real-time  inputs  (10  points  per  second)  to  the  SURC  equations  where  the 
program  was  set  to  run  with  no  layer  present.  (See  program  description, 
Appendix  III) . 

Figure  6  shows  the  difference  between  the  corrected  elevation  angle 
data  where  the  AFWTR  data  in  all  cases  is  subtracted  from  the  SURC  data. 
During  the  first  10  seconds  of  flight  the  true  elevation  angle  difference  (B0) 
reaches  a  peak  of  about  0.  0075  degrees  (0.  131  mr).  This  difference  decreases 
rapidly  to  a  negligible  amount  after  the  first  50  seconds  of  flight.  Figure  7 
shows  the  elevation  angle  error  data,  €,  obtained  with  the  AFWTR  ray-tracing 
program  and  figure  8,  the  same  angle  error  data  from  the  SURC  equations. 

It  is  apparent  that  the  largest  angle  difference  (0.0075  degrees)  in  figure  6  is 
contributed  by  a  perturbation  in  the  AFWTR  calculations.  The  refractivity 
profile,  figure  5,  does  not  show  any  discontinuity  above  the  site  elevation 
and  the  vehicles  elevation  and  range  are  mono  tonically  increasing  (figures  3 
and  4).  Therefore,  the  perturbation  shown  in  the  AFWTR  data  appears  to  be 
due  to  some  computational  error  in  the  ray-tracing  program  and  not  pro¬ 
duced  because  of  the  propagation  conditions. 

Figure  9  shows  the  differences  between  the  corrected  ranges.  These 
differences  are  so  small  (well  below  the  range  resolution  of  the  radar)  they 
can  be  considered  to  be  negligible.  These  results  show  that  the  exponential 
model  atmosphere  does  provide  an  accurate  description  of  the  radio  propa¬ 
gation  conditions  for  this  particular  set  of  tracking  conditions. 
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FIGURE  4.  RANGE  VS.  TIME 
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FIGURE  5.  RADIO  REFRACTIVITY  PROFILE  (SOUTH  VANDENBERG) 
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FIGURE  7.  ELEVATION  ANGLE  ERROR  VS. TIME  (AFWTR) 
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FIGURE  8.  ELEVATION  ANGLE  ERROR  VS.  TIME  (SURC) 
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Figures  10  and  11  show  the  tracking  parameters  from  the  radar 
(003003)  at  Point  Mugu.  Figure  12  shows  the  difference  between  the  SURC 
and  AFWTR  elevation  angle  data.  Initially,  for  tracking  angles  below  five 
degrees  (figure  10),  the  difference  is  about  .015  degrees  (0.26  mr)  and 
thereafter  decreases  below  0.  003  degrees  (0.  053  mr).  Figure  13  shows  the 
corrected  range  error  difference  and  the  maximum  difference  reaches  about 
seven  feet  at  low  tracking  angles  and  then  decreases  below  1.  5  feet.  Target 
position  is  most  seriously  affected  by  the  angle  errors  and  the  angle  dif¬ 
ference  (figure  12)  of  0.  015  degrees  at  a  range  of  500  KFT  represents  a 
height  difference  of  130  feet.  At  this  time,  the  vehicle  is  around  30  KFT 
altitude.  At  200  seconds  the  height  difference  is  only  98  feet  at  a  vehicle 
altitude  of  about  one  million  feet. 

Figures  14  through  17  show  corresponding  data  for  the  Point  Pillar 
radar  (2  13001).  This  radar  picks  the  vehicle  up  at  a  higher  initial  tracking 
angle  and  range  due  to  the  fact  it  is  located  far  north  of  the  firing  site.  The 
angle  and  range  difference  data  is  obviously  negligible. 

3 .  1  Some  Comments  on  the  Analysis  of  Operation  5265 

Except  for  tracking  elevation  angles  below  five  degrees  the  SURC 
real-time  calculations  are  essentially  in  agreement  with  the  post-flight,  ray- 
traced,  results  provided  by  AFWTR.  Of  course,  the  greater  differences  in 
the  elevation  angle  comparisons  which  occur  at  short  ranges  do  not  produce 
large  position  errors  in  height.  These  angle  differences  would  be  serious  at 
long  ranges  if  they  remained  above  0.006  degrees  (0.  1  mr).  A  more  interest¬ 
ing  comparison  could  be  made  when  the  radars  are  tracking  at  small  eleva¬ 
tion  angles  and  very  long  ranges  as  the  vehicle  goes  over  the  radar  horizon. 
Unfortunately,  the  tracking  data  for  operation  5265  terminated  long  before 
the  radar  horizon  was  reached. 

A  more  interesting  comparison  is  made  in  the  analysis  of  Operation 
2062,  which  follows,  where  one  radar  does  track  the  vehicle  almost  to  the 
radar  horizon. 
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FIGURE  10.  ELEVATION  ANGLE  VS.  TIME  (OP5265  PT.  MUGU) 
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FIGURE  12.  DIFFERENTIAL  CORRECTED  ELEVATION  ANGLES  VS.  TIME 
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FIGURE  14.  ELEVATION  ANGLE  VS.  TIME  (OP5265  PT.  PILLAR) 
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SECTION  IV 


COMPARISON  OF  TRACKING  ERRORS  BETWEEN  SURC 
CALCULATIONS  AND  AFWTR  FAY-TRACED  RESULTS  - 

OPERATION  2062 

Figures  18  through  23  show  the  tracking  and  difference  data  for 
Operation  2062  for  the  South  Vandenberg  tracker  (023002).  From  figure  21 
it  may  be  noted  that  the  corrected  elevation  angle  data  obtained  from  the 
AFWTR  ray-traced  program  is  perturbed  within  the  first  few  seconds  of 
flight  Figure  24  shows  the  associated  radiosonde  profile  and  it  is  apparent 
that  the  propagation  conditions  do  not  suggest  a  reason  for  this  behavior.  As 
noted  before,  the  AFWTR  calculations  of  the  elevation  angle  error  (figure  21) 
during  the  first  few  seconds  of  flight  are  not  compatible  with  the  smooth  in¬ 
crease  of  elevation  angle  and  range  (figures  18  and  19). 

The  difference  between  SURC  and  AFWTR  calculations  of  the  cor¬ 
rected  elevation  angle  (figure  22)  and  range  (figure  23)  are  negligible  after 
these  first  few  seconds  of  flight. 

Figures  25  through  28  show  equivalent  graphical  results  with  the 
tracking  radar  located  at  Point  Mugu  (003001).  As  before,  the  difference 
between  corrected  elevation  angles  and  ranges  (figures  27  and  28)  are 
negligible  beyond  30  seconds  of  the  flight.  The  radiosonde  data  (Appendix  II) 
indicates  a  very  weak  layer  might  be  present  around  3  KFT.  Since  the  layer 
was  not  introduced  in  the  SURC  equations  one  might  expect  that  the  elevation 
angle  error,  would  be  larger  using  the  AFWTR  program.  Therefore, 
during  the  first  few  seconds  of  tracking  from  Point  Mugu  the  effect  of  a  layer 
would  cause  the  elevation  angle  difference  (figure  27)  to  be  larger  The 
results  shown  are,  therefore,  consistent  with  such  propagation  conditions. 

Figures  29  through  32  show  the  results  obtained  with  the  tracker 
located  at  Point  Pillar  (213001).  The  comparative  accuracies  obtained 
throughout  the  first  300  seconds  of  tracking  results  are  very  good.  The 
initial  range  error  difference  of  about  six  feet  (figure  23)  is  slightly  larger 
than  in  previous  cases.  However,  it  should  be  noted  that  this  time  the 
vehicle  was  at  a  range  of  about  one  million  feet  (approximately  200  miles). 

Figures  33  through  36  show  results  for  the  next  300  seconds  of  flight 
from  Point  Pillar  where  the  vehicle  is  now  heading  towards  the  radar  horizon. 
From  figure  33  it  may  be  seen  that  the  elevation  angle  is  decreasing  towards 
five  degrees  and  the  range  (figure  34)  increasing  beyond  10  million  feet 
(approximately  2000  miles).  A  loss  of  signal  occurs  momentarily  around 
530  seconds. 
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FIGURE  18.  ELEVATION  ANGLE  VS.  TIME  (OP2062  SVBG) 
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FIGURE  19.  RANGE  VS.  TIME 
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FIGURE  20.  ELEVATION  ANGLE  ERROR  VS.  TIME  (SURC) 
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FIGURE  21.  ELEVATION  ANGLE  ERROR  VS.  TIME  (AFWTR) 
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FIGURE  22.  DIFFERENTIAL  CORRECTED  ELEVATION  ANGLES  VS.  TIME 
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Figure  24.  RADIO  REFRACTIVITY  PROFILE  (SOUTH  VANDENBERG) 
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FIGURE  25.  ELEVAliUJiN  ANGLE  VS.  TIME  (OP2062  PT.  MUGU) 
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FIGURE  26.  RANGE  VS.  TIME 
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FIGURE  29.  ELEVATION  ANGLE  VS.  TIME  (OP2062  PT.  PILLAR) 
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FIGURE  30.  RANGE  VS.  TIME 
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FIGURE  32.  DIFFERENTIAL  CORRECTED  RANGES  VS.  TIME 
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As  the  tracking  elevation  angle  decreases  towards  five  degrees  the 
elevation  angle  difference  (figure  35)  increases  to  a  value  around  0.003 
degrees  (0.  0525  mr)  which  is  still  less  than  the  reported  angular  resolution 
of  the  tracker.  The  range  difference  (figure  36)  increases  towards  a  value 
of  about  five  feet  at  a  vehicle  range  of  1  0  million  feet. 

4 .  1  Some  Comments  on  the  Analysis  of  Operation  2062 

Basically,  the  analytical  results  of  Operation  2 062  confirm  that  the 
SUR.C  equations  can  provide  refraction-induced  error  data  in  agreement  with 
ray-traced  results.  As  the  tracking  elevation  angle  decreases  below  five 
degrees  the  differences  between  the  radiosonde  refractivity  measurements 
and  the  CRPL  exponential  atmosphere  become  more  apparent.  This  subject 
will  be  treated  in  greater  detail  in  the  conclusions. 

It  is  noteworthy  that  the  agreement  between  methods  is  also  very  good 
at  long  tracking  ranges  as  shown  by  the  final  set  of  results.  As  the  elevation 
angle  from  the  Point  Pillar  tracker  decreased  below  five  degrees  the  tracking 
data  became  erratic  probably  due  to  multipath  effects.  Therefore,  the  plots 
were  terminated  at  600  seconds. 
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SECTION  V 


COMPARISON  OF  TRACKING  ERRORS  BETWEEN  SURC 

CALCULATIONS  AND  GERTS  FLIGHT  DATA  -  OPERATION  7335 

Figures  37  and  38  show  the  elevation  angle  and  range  of  the  vehicle 
from  the  GERTS  system  ( Vandenberg) .  Figure  39  shows  the  elevation  angle 
error  calculated  with  the  SURC  equations  where  from  Appendix  II  the  as¬ 
sociated  radiosonde  data  did  not  indicate  any  elevated  layer  or  non-standard 
ref ract  ivity  behavior  However,  figure  40,  which  shows  the  GERTS  elevation 
angle  error  calculations,  indicates  a  peak  in  the  region  of  5  seconds.  This 
behavior  is  difficult  to  explain  since  both  range  and  elevation  angle  are 
monotonically  increasing 

Figure  41  shows  the  difference  between  the  true  elevation  angle  cal¬ 
culations  and  the  minimum  around  5  seconds  is  clearly  due  to  the  GERTS 
calculation  described  above.  The  angle  difference  decreases  rapidly  to  a 
negligible  amount  after  25  seconds  The  initial  maximum  difference  of  about 
-0.  010  degrees  is  comparable  with  the  results  of  the  previous  operation. 

Figure  42  shows  the  corrected  range  difference  which  decreases 
rapidly  after  the  first  20  seconds  of  flight  Since  the  GERTS  tracker  is 
particularly  interested  in  range-rate  data  it  is  significant  to  note  that  the 
difference  in  calculation  is  not  only  small  but  remains  essentially  constant 
after  forty  seconds  of  flight. 
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FIGURE  38.  RANGE  VS.  TIME 
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FIGURE  39.  ELEVATION  ANGLE  ERROR  VS.TIME  (SURC) 
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FIGURE  40.  ELEVATION  ANGLE  ERROR  VS.  TIME  (GERTS) 
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SECTION  VI 


SOME  COMMENTS  ON  THE  GERTS  REAL-TIME  CORRECTION 

METHOD 

The  GERTS  radar  data  was  corrected  using  a  real-time  method. 

These  equations  are  derived  from  Snell’s  law  for  refraction  in  spherical 
coordinates.  4  To  calculate  the  true  elevation  angle  an  expression  is  used, 
where 

n  sin  E0'=  n0  sin  E0.  (12) 

and 

En'  =  the  true  elevation  angle 

=  the  measured  elevation  angle 
n0  =  the  surface  refractive  index 

n  =  the  average  index  of  refraction  between  the  missile  and  the 
tracker . 

This  equation  can  be  developed  using  small  angle  approximations  to 
give  an  expression  for  the  elevation  angle  error,  AE,  in  microradians, where 

£  E  =  tan  E0  (N0  -  N)  (13) 

and 

N0  =  the  surface  refraetivity 
=  (  n,  -  1)  x  10* 

and 

N  =  the  average  index  of  refraetivity  between  the  missile  and  the 
radar,  where 

N  =  (n  -  1 )  x  1 06 

The  problem  is  to  find  values  for  N  which  provide  meaningful  values 
for  the  refraction-induced  error  data.  By  assuming  an  exponential  re- 
fractivity  model  of  the  CRPL  type2  the  average  refraetivity  over  the  path  is 
approximated  by 


where  Cp  is. the  surface  dependent  exponential  scale  height  and  the  target 
height  H  is  found  from  apparent  radar  data. 
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There  are  certain  limitations  imposed  by  this  method  of  development 
which  can  limit  the  usefulness  of  these  real-time  equations. 

In  the  first  place,  Snell's  Law  is  not  expressed  in  terms  of  the  average 
of  refractivity  (12), but  it  defines  the  relationship  of  the  index  of  refraction 
and  the  local  elevation  angle  of  a  ray  at  a  point  in  space.  Second,  the  target 
height  used  in  (14)  is  calculated  from  apparent  radar  data  which  is  meaning¬ 
ful  only  at  large  elevation  angles.  The  calculated  elevation  angle  errors 
using  these  equations  must,  therefore,  be  in  error.  By  relating  the  data 
with  accurate  ray-traced  results  certain  empirical  corrections  can  be  added 
to  these  equations  to  obtain  meaningful  data,  which  is  apparently  what  has 
been  done. 

The  comparison  between  these  calculations  and  the  SURC  results, 
(figure  41),  indicates  that  the  GERTS  equations  provide  very  good  error  data 
for  this  particular  operation. 

The  major  difficulty  will  be  realized  if  it  is  required  to  extend  the  use 
of  these  equations  to  low  elevation  angles  and  to  provide  meaningful  data 
when  a  strong  inversion  is  present. 

The  most  serious  limitation  is  that  the  doppler  error  angle,  6, 

(figure  1)  cannot  be  calculated  from  the  GERTS  equation,  where, 

6  =  r  -  c.  (15) 

These  equations  do  not  permit  a  calculation  of  the  ray  path  bending,  T.  Since 
the  doppler  error  angle,  6,  corrects  velocity  measurements  obtained  from 
doppler  or  range-rate  measurements  a  knowledge  of  this  angle  is  necessary 
to  obtain  accurate  tangential  velocity  information  at  low  tracking  angles.3 
This  situation  occurs  when  a  high  trajectory  object  is  being  tracked  at  low 
elevation  angles  over  the  radar  horizon. 

It  is  sufficient  to  note  that  the  SURC  equations  have  none  of  these 
limitations  since  they  are  derived  in  direct  relationship  to  the  physical 
principles  of  ray  behavior  in  a  refractive  medium,  including  the  effect  of  an 
elevated  layer.  It  is,  therefore,  advantageous  to  use  these  SURC  equations 
if  their  execution  time  can  be  made  compatible  with  the  overall  mission 
requirements.  These  equations  can  be  directly  extended  for  use  in  the 
optical  tracking  region  should  such  a  requirement  arise. 
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SECTION  VII 


SUMMARY  COMMENTS 

The  analysis  of  three  AFWTR  missions  demonstrates  that  the  SURC 
equations  can  be  used  in  near  -  real-time  to  correct  ref  raction  -  indue  ed  tracking 
errors.  The  corrected  tracking  data  using  these  equations  is  comparable  in 
accuracy  with  the  use  of  ray-tracing  correction  methods.  They  can  be  used 
with  a  minimum  of  meteorological  support  and  can  take  account  of  additional 
effects  produced  by  an  elevated  layer.  This  latter  ability  gives  them  a 
distinct  advantage  over  the  present  GERTS  real-time  correction  method. 

Also,  the  GERTS  equations  do  not  provide  a  correction  to  velocity  errors, 
which  are  introduced  because  of  ray -bending.  This  inability  may  become 
significant  if  greater  velocity  accuracy  is  required  in  future  missions. 

Since  any  computational  method  requires  some  finite  execution  time, 
there  is  no  direct  method  to  obtain  truly  real-time  error  data.  However, 
the  SURC  equations  can  be  programmed  for  use  with  a  raw  radar  data  rate  of 
greater  than  ten  points  per  second.  A  simple  linear  extrapolation  of  the 
error  data  would  permit  equivalent  real-time  error  data  to  be  made  available. 
Since  the  required  extrapolation  period  is  set  by  the  particular  processing 
machine  we  have  not  attempted  to  introduce  this  relatively  simple  function 
in  the  present  program. 

The  question  of  the  magnitude  of  residual  tracking  errors  following 
corrections  for  refraction  was  covered  in  the  previous  report.1  Unless  the 
two  dimensional  structure  of  refractivity  in  height  and  distance  could  be 
measured  instantaneously,  the  re  is  no  analytical  method  to  determine  the 
magnitude  of  these  residual  errors.  For  example,  figure  43  shows  the 
nature  of  the  horizontal  changes  in  refractivity  which  were  measured  with  a 
refractometer-equipped  aircraft.  The  vertical  variations  of  refractivity 
near  the  antenna  have  the  greatest  effect  on  the  magnitude  of  height  errors. 
However,  variations  progressively  further  away  from  the  radar  contribute  in 
a  range -weighted  manner  to  the  final  overall  height  errors.  Therefore,  in¬ 
creasing  the  detail  of  refractivity  structure  near  the  radar  site  will  not 
necessarily  decrease  the  overall  height  error  directly 

A  sensible  approach  to  the  determination  of  residual  errors  in  real¬ 
time  tracking  operations  is  to  confront  the  corrected  radar  position  data  with 
independent  measurements  of  target  position.  As  mentioned  in  the  introduction, 
radar-altimeter  equipped  aircraft  and  established  satellite  orbital  parameters 
could  be  used  to  fix  the  target  locations. 
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Measurements  of  this  kind  have  already  been  carried  out  with  aircraft6 
where  a  sea- reflection  interferometer  was  used  to  measure  the  angle  of 
arrival  of  the  signal.  In  general,  the  receiving  system  must  be  able  to  re¬ 
solve  angles  of  arrival  of  at  least  0.  05  milliradians  to  be  able  to  discriminate 
the  effect  of  residual  refraction  errors.  This  accuracy  requirement  essen¬ 
tially  rules  out  using  the  FPS-16  radars  in  such  an  experiment.6 

The  exponential  model  atmosphere  including  the  effect  of  the  in¬ 
version  layer  has  demonstrated  an  ability  to  represent  propagation  conditions 
by  comparison  with  ray-traced  results.  The  effect  of  the  layer  is  covered  in 
Appendix  I  .  It  appears  from  analysis  covered  in  the  previous  report1  that  in 
the  climatological  environment  peculiar  to  the  coast  of  California,  that  the 
CRPL  exponential  could  be  adjusted  to  provide  better  agreement  with  ray- 
traced  results.  By  establishing  the  relationship  between  surface  values  and 
the  initial  gradient  over  the  first  one  kilometer  height  section,  a  regression 
analysis  can  be  carried  out  using  available  radiosonde  data.  The  indication 
based  on  previous  analysis  is  that  the  exponential  decay  constant,  c,  used 
with  the  exponential  model  should  be  slightly  smaller  for  the  California 
coastal  area  than  is  given  by  the  CRPL  model.  To  this  end,  data  which  has 
been  collected  by  ray-tracings  of  radiosonde  and  aircraft  measurements  is 
compiled  in  an  addendum  of  which  a  limited  number  of  copies  are  available 
for  this  type  of  further  analysis.  Since  this  type  of  information  is  not  of 
general  interest,  the  data  is  not  included  with  this  report. 
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SECTION  VIII 


REC  OMMENDATIONS 

1.  Since  the  SURC  refraction-correction  program  can  be  used  in  near- 
real-time  to  correct  range,  angle  and  velocity  data,  it  is  recommended  that 
it  be  initially  implemented  on  the  AFWTR  SDS-930  computer  for  further 
evaluation.  The  program  was  run  on  the  SURC  SDS-930  using  AFWTR  radar 
data  tapes.  Additional  evaluation  of  the  program  should  be  made  by  AFWTR 
particularly  when  tracking  is  performed  in  the  presence  of  a  strong  elevated 
layer  (inversion). 

2.  The  program  operators  should  become  completely  confident  with  the 
program  and  evaluate  it  under  various  optional  methods  of  inputing  meteoro¬ 
logical  data.  It  is  then  recommended  that  an  evaluation  be  made  of  the 
reduced  running  time  possible  through  the  use  of  machine  language  operation 
on  other  superior  software. 

3.  An  evaluation  of  the  limitations  of  the  electromechanical  angle  accuracy 
of  the  FPS-16  trackers  should  be  carried  out  to  determine  if  they  can  be  used 
to  experimentally  determine  the  residual  errors  after  real-time  corrections 
are  applied.  If  the  accuracy  is  well  under  0.  1  milliradians  it  would  be  use¬ 
ful  to  compare  the  corrected  radar  data  using  the  SURC  program  with 
independently  determined  target  position  data.  Airborne  radar  altimeters 
should  be  able  to  give  target  (aircraft)  heights  to  within  50  feet.  The  pos¬ 
sibility  of  using  satellite  targets  should  be  considered. 

4.  Presently  available  radiosonde  data  can  be  used  to  perform  a  cor¬ 
relation  analysis  of  the  initial  1  km  gradient  with  surface  data.  In  this  way 
it  appears  possible  to  improve  the  correction  available  with  the  SURC  pro¬ 
gram  by  determining  an  exponential  model  more  suited  to  the  seasonal 
characteristics  in  the  Vandenberg  area. 

5.  It  is  recommended  that  AFWTR  consider  using  the  consultation  of 
SURC  in  these  areas  in  view  of  our  present  experience  with  the  nature  of 
refraction  errors  affecting  the  AFWTR  tracking  system. 
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APPENDIX  I 


A  METHOD  TO  INCLUDE  THE  EFFECT  OF  AN  ELEVATED  LAYER 

AND  TEST  FOR  TRAPPING 


1.0  A  METHOD  TO  INCLUDE  THE  EFFECT  OF  AN  ELEVATED  LAYER 
AND  TEST  FOR  TRAPPING 


On  the  average  the  exponential  model  atmosphere  represents  realistic 
conditions  because  of  the  exponential  decrease  of  pressure  with  height.  The 
model  becomes  very  non-standard  when  an  inversion  exists,  in  which  case 
the  water  vapor  tends  to  remain  under  the  inversion  and  then  rapidly  de¬ 
creases  with  height  through  the  inversion  layer.  In  terms  of  refr ac tivity , 

N,  is 


N  L 


77.  6 
T 


(P  + 


48 1  Oe 
T  ) 


(1A) 


where 


p  =  the  air  pressure  (mb) 

T  =  the  air  temperature  (°  K) 

e  =  the  associated  water  vapor  pressure  (mb) 

The  rapid  decrease  of  water  vapor  through  the  inversion  causes  N  to 
also  decrease  thereby  forming  an  elevated  layer.  This  layer  causes  the 
ray-bending  to  increase  above  normal  when  the  ray  enters  the  layer  at  small 
angles  of  incidence.  If  the  layer  produces  enough  bending  the  ray  does  not 
immediately  pass  through  the  layer  and  a  radio  hole  is  formed. 

Figure  1-1  shows  two  refractivity  profiles,  one  the  CRPL  exponential 
model  and  the  other  typical  of  propagation  conditions  in  the  presence  of  a 
strong  inversion.  The  CRPL  exponential  is  generated  using  the  surface 
value  of  345  N  units. 


Figure  1-2  depicts  the  paths  of  two  rays,  one,  S0,  propagated  in  an 
exponential  atmosphere  and  the  other,  Sj^,  in  the  same  atmosphere  but  with 
an  elevated  layer  included.  Due  to  the  increased  refractivity  gradient  in  the 
inversion  layer  the  ray,  S^,  is  bent  significantly  as  it  passes  through  the  layer. 


The  bending  in  the  layer  produces  a  local  elevation  angle,  9P,  at  the 
top  of  the  layer,  where 


9P-  [0T  + 


2DH 


2DN  -.1 
1 0B  J 


(2A) 


(see  Radio  Meteorology,  page  83,  NBS  Monograph  92), 
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HEIGHT  (km) 


60 


FIGURE  1-1.  THE  EXPONENTIAL  PROFILE  AND  A  LAYERED  MODEL 
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FIGURE  1-2.  RAY  TRACING  GEOMETRY  IN  A  LAYERED  MODEL 
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whe  r  e 


=  the  local  angle  of  incidence  of  the  ray  into  the  layer 

DH  =  the  layer  thickness  (positive) 

DN  =  the  change  in  N  through  the  layer,  and  is  defined  here  as  positive. 

Since  the  inversion  tends  to  trap  the  water  vapor  under  the  layer  the 
concentration  of  water  vapor  with  height  decreases  slowly  under  the  layer. 
Examination  of  meteorological  data  shows  that  the  ref ractivity ,  N,  then 
tends  to  decrease  almost  linearly  with  height  under  a  strong  inversion.  It  is, 
therefore,  convenient  to  obtain  a  determination  of  9V  and  the  distance*  RC,  to  the 
layer,  RC ,  by  using  an  effective  earth  radius  approach. 


Referring  to  figure  1-2,  the  true  earth's  radius  is  replaced  by  an 
effective  earth  radius,  AE,  where  (Schelleng,  et  al,  Proc,  IRE  21 ,  March  1933). 


AE 


1  + 


I0e 


-  Hs 


The  angle  of  incidence,  0T,  is  then  given  by 
0T  =  tt/2  -  rv 

whe  re 


(3A) 


(4A) 


ry 


sin"1  [cos  9n  • 


AE  +  h: s , 
AE  +  HI  J 


The  range,  RC,  where  the  ray  enters  the  layer  is  then 
RC  =  (AE  +  HI)  •  cos  rv  .  (AE  +  HS)  sin  90. 


(5A) 


(6A) 


It  is  mathematically  convenient  to  concentrate  the  bending  produced 
by  the  layer  on  Sl  at  the  point  where  the  ray  enters  the  layer.  In  other 
words,  as  the  ray  Sj^  enters  the  layer  at  an  angle,  9 T,  assume  that  it 
abruptly  changes  its  direction  to  a  new  angle,  9p,  given  by  (2A).  Now  com¬ 
pare  the  remaining  bending,  T,  between  Se  and  S-^  over  the  remaining  range 
R  -  RC.  If  the  value  of  ref  ractivity,  N^,  at  the  target  is  very  small  com¬ 
pared  to  the  value,  Nv  at  the  base  of  the  layer,  then,* 

T  T  Nn  cot  p.  X  10"6  (7A) 

where  P.  is  the  initial  elevation  angle  of  the  ray  at  the  base  of  the  layer. 

*Bean,  B.R.  and  E.  J.  Dutton,  "Radio  Meteorology,  M  NBS  Monograph  92, 
March  1966,  p.  53. 
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(8A) 


Therefore,  over  the  remaining  path,  R  -  RC, 

T l  •  cot  0^  x  10“6  for  S-^ 

where  N1  is  the  value  of  refractivity  at  the  base  of  the  layer  obtained  from 
radiosonde  data  (figure  I-l)  and  for  the  other  ray  Se  in  the  exponential 
atmosphere, 

Te  -N1  •  cot  ei  X  10'6  (9A) 


Taking  the  ratio  of  (8A)  and  (9A)  the  bending  of  the  ray  over  the  remainder 
of  the  path,  R  -  RC,  is  changed  from  its  normal  bending,  Tp>  in  the  ex¬ 
ponential  atmosphere  where 


~  cotj^ 

rL  Te  cot 


(10A) 


Since  the  effect  of  bending  is  only  significant  at  small  angles  of  incidence 
then  using  small  angle  approximations 


L  ci 


0 


( 1 1  A) 


where  Te  is  obtained  directly  from  (3)  over  the  path  where  the  range,  R, 
exceeds,  RC . 


Therefore, 

.  (8) 

(Section  2.  0),  can  be  redefined  when  a  layer  is 

present  by  setting 

n 

n 

T 

sin  (fl0 

-  Ti) 

6R^  =  23  sin  (0O  -  T[  •  m)  6R{ 

(12A) 

i=l 

i=  1 

where 

m 

=  1  for  R  <RC 

(13A) 

m 

II 

CD 

<D 

for  R  >  RC 

(14A) 
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2.0 


A  COMPARISON  OF  CALCULATIONS  INCLUDING  A  LAYER  WITH 
RAY -TRACING  ANALYSIS 


Let  the  radiosonde  data  of  figure  .1-1  describe  propagation  conditions 
in  the  presence  of  a  strong  inversion.  Figure  1-3  shows  a  comparison  of 
results  with  the  layer  present  between  the  SURC  calculations  and  ray-traced 
results  where  the  tracking  station  height  is  at  0.  5  km  above  sea  level.  With 
the  layer  not  included  the  error  angle  is  seen  to  be  significantly  less  at  small 
elevation  angles. 

The  elevation  angle  error,  €,  from  SURC  calculations  is  given  by 

i  n 

f  =  %  -  sin_1  [r  S  sin  (qo  -  Tj  •  m)  6rJ  (15A) 

i=  1 

whe  re 

m  =  1  for  R  <  RC  (l6A) 

m  =  for  R  >RC  (17A) 

and  RC,  and  %  are  determined  from  (6A),  (4A),  and  (2A),  respectively, 
and  T ^  by  (3),  The  range  intervals,  6Rj,  were  found  from  (10)  using  nine 
intervals  or  N  equals  10  in  equation  (11) 

For  the  station  height  of  0.5  km  the  station  value,  Ns  is  313  and  all 
other  parameters,  c,  k0,  kx,  kP,  kQ,  are  related  to  this  value  of  Ns  as  pre¬ 
viously  reported. 

Figure  I  -4  shows  a  comparison  where  the  station  height  is  increased 
to  1  km  above  the  surface.  As  the  distance  between  the  base  of  the  inver¬ 
sion  and  the  station  dec  reases,  the  effect  of  the  layer  on  ray -bending  be¬ 
comes  greater.  The  calculations  with  the  layer  excluded  (that  is,  m  is 

set  equal  to  unity  for  all  ranges  (16A))  shows  very  much  less  angle  error  at 
small  elevation  angles.  This  data  is  equivalent  to  an  exponential  atmosphere. 

Obviously,  including  the  effect  of  the  layer  in  the  SURC  equations 
gives  a  first  order  improvement  to  the  prediction  of  bending  errors  when  a 
layer  is  present.  It  is  also  apparent  and  previously  reported  that  even  the 
effect  of  a  strong  layer  has  relatively  little  effect  on  ray -bending  until  the 
elevation  angle  decreases  below  two  degrees  (35  mr). 
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FIGURE  1-4.  A  COMPARISON  OF  ELEVATION  ANGLE  ERRORS 


Referring  to  (10)  the  range  increments,  6R^,  are  fairly  large  in  which 
case  the  effect  of  the  layer  is  generally  not  introduced  into  (15A)  until  after 
the  range,  R,  exceeds  RC .  This  effect  compensates  for  the  assumption  made 
in  the  derivation  that  the  bending  produced  by  the  layer  is  concentrated  at 
the  point  where  the  ray  enters  the  layer.  Also,  the  assumption  that  the  ray¬ 
bending  is  directly  proportional  to,  cot  6.,  (7A)  is  true  only  when  the  range 
is  very  large.  Therefore,  the  above  comparative  results  (figures  I -3  and 
1-4)  are  most  applicable  when  tracking  a  vehicle  near  the  radar  horizon. 
However,  the  target  height  error,  for  a  given  elevation  angle  error,  €,  is 
almost  directly  proportional  to  range.  That  is,  the  height  error,  Ah,  is 

Ah  ~  f  •  R  ( 18A) 

Therefore,  the  errors  in  calculating  the  effect  of  a  layer  from  (15A)  at  short 
tracking  ranges  is  much  less  serious  than  if  these  errors  were  implicit  in 
the  calculations  for  long  ranges. 


To  increase  the  accuracy  of  these  calculations  at  short  range  would 
require  using  the  expression  for  from  (3),  then 

—  Na  cot  v2  x  10"6  [l  -  — ^ -  e  kl  Zkigi]  (19A) 

m  +  Si 

where  ga,  and  kl  are  associated  with  the  angle,  02 . 

Also 

Tp  -  N  cot  VT  x  10"6  [1  -  - - •  e"ks2  '  2k?gs]  (20A) 

e  1  1  ks  +  gP 


where  Vx ,  g?,  and  k?  are  associated  with  the  angle  9^.  S  ince  the  range, 

R  -  RC,  is  equal  in  (19A)  and  (20A),  then  neglecting  second  order  effects, 


TL  “  T 


e 


cot  yo 
cot  Vx 


-  T, 


Vj 


(21A) 


For  example,  Vx,  for  small  elevation  angles  is  essentially 


Va  “  0X  +  (k0  +  ka  e 


-k3  (R  -  RC) 


( 22A) 
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For  real-time  application  where  these  SURC  equations  are  most 
useful,  complicating  the  summation  in  (15A)  through  the  use  of  equation 
(21A)  is  not  justified. 

Furthermore,  the  characteristics  of  the  layer  generally  change 
significantly  with  distance  from  the  radar.  The  radiosonde  measurement 
of  the  layer  characteristics  near  the  site  may  be  unrelated  to  the  layer 
conditions  at  the  range  RC.  Therefore,  including  the  layer  in  the  SURC 
equations  only  provides  a  first-order  correction  in  any  case  and  in  the  same 
way  even  if  ray-tracing  calculations  are  used  the  same  problems  in  defining 
the  layer  apply.  However,  for  very  small  elevation  angles  the  elevation 
angle  error,  c,  can  be  partly  compensated  by  including  the  factor,  01/0?f  in 
the  bending  equation  when  the  range  exceeds,  RC, 
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A  TEST  FOR  TRAPPING  CONDITIONS 


3.  0 


When  the  measured  characteristics  of  a  layer  are  introduced  the  pre¬ 
ceding  analysis  shows  th&t  a  first-order  correction  to  ray-bending  can  be 
made.  Possibly  of  equal  importance  is  the  fact  that  these  equations  can  be 
used  to  give  a  prediction  of  the  tracking  elevation  angle  where  trapping  could 
occur.  For  example,  setting  equal  to  zero  in  (2A)  defines  a  limiting 
condition  where  the  ray  would  be  trapped  for  a  value  of  less  than 


2DN 

106‘ 


2DH 


(23A) 


Obviously,  since  0T  has  to  be  a  real  number  the  first  condition  to  produce 
trapping  is  that 


DN  ^  DH 

i3n- 


(24A) 


If  these  conditions  are  satisfied  then  0O  is  found  from  an  expanded  form  of 
Snell's  law  in  polar  coordinates  where 

ns  rg  cos  0Q  =  n1  rL  cos  0T  (25A) 

and  since 

n  s  =  If"  Ns/106 

rs  =  r0  +  Hs 

nl  =  1  +  N1/ 1 06 


neglecting  second  order  effects 


cos  0orp  — 


COS  01 


h  +  Ns_-  Ni 

L1  +  106 


H,  -  Hc 


(26A) 


For  a  real  value  00^  the  right  hand  side  of  (26A)  must  be  less  than  one.  If 
it  i^  then  0O^  is  determined  and  represents  the  largest  apparent  elevation 
angle  which  could  be  trapped.  The  radar  operator  is  then  alerted  to  the 
fact  that  the  target  signal  may  be  lost  for  tracking  angles  less  than  0O^. 
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With  the  radar  at  0.  5  km  above  sea  level  both  the  ray-tracing 
analysis  and  the  SURC  equation  (26A)  indicate  that  the  trapping  angle, 
was  zero  degrees,  or  in  other  words,  no  ray  would  be  trapped  for  apparent 
elevation  angles  above  zero  degrees.  With  the  radar  at  1  km  above  sea  level 
ray-tracing  results  show,0orp  equals  0.2817  degrees  (4.92  mr)  and  the 
SURC  equations  give  0.2480  degrees  (4.335  mr).  Such  close  agreement 
indicates  that  the  SURC  equations  can  be  used  to  give  a  very  good  estimate 
of  the  trapping  angle,  0o,p,  compared  to  ray-tracing  methods.  However, 
again  it  should  be  emphasized  that  whether  or  not  these  results  are  meaning¬ 
ful  depends  upon  the  spatial  characteristics  of  the  layer. 
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APPENDIX  II 
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APPENDIX  III 


DESCRIPTION  OF  THE  SURC  REAL-TIME 
REFRACTION  CORRECTION  TECHNIQUE 

1. 0  INTRODUCTION 

The  enclosed  program  used  on  the  SDS-930  computer  at  SURC  con¬ 
sists  of  several  sections,  each  section  individually  numbered.  For  ease  in 
explanation,  the  total  package  is  also  renumbered  by  consecutive  pages 
from  94  through  120. 

It  should  be  pointed  out  that  the  actual  real-time  refraction  correc¬ 
tion  program  is  described  in  pages  109  through  112  which  is  a  small  part  of 
the  program  described  herein. 

The  total  package  is  actually  two  individual  programs: 

1.  The  initialization  program,  pages  94  through  108,  reads  the 
data  cards,  calculates  some  a  priori  parameters  and  outputs 
them  on  magnetic  tape. 

2.  The  actual  refraction  correction  program,  pages  109  through 
120,  reads  in  a  priori  parameters  and  the  radar  data,  cal¬ 
culates  the  refraction  error  and  outputs  the  results  on  either 
the  printer  or  magnetic  tape,  whichever  is  selected  by  sense 
switch  control. 
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2.0 


THE  INITIALIZATION  PROGRAM 


2 .  1  Pages  94  -  96  :  Refraction  Correction  Initialization 

This  program  computes  coefficients  which  are  required  for  use  in  the 
real-time  running  program  (Refraction  Correction  -  page  109).  Meteoro¬ 
logical  information  is  entered  at  this  point  as  well  as  the  height  of  the  radar 
above  mean  sea  level,  HS,  in  kilomete  rs . 

Several  options  are  available  for  entry  of  meteorological  data: 


If  the  radar  station  value  of  refractivity,  SNS,  is  known,  enter 
it  directly.  In  general  this  value  would  be  obtained  from  a 
radiosonde  profile  where  the  value  SNS  is  picked  from  the 
profile  at  the  height,  HS,  of  the  particular  radar  above  mean 
sea  level.  The  radiosonde  profile  would  be  selected  from  the 
local  launch  station  and  as  close  in  time  to  the  launch  as 
possible. 

Using  this  method  the  radiosonde  profile  will  indicate  the 
presence  or  absence  of  the  inversion  layer.  If  the  layer  is 
present  and  above  the  height  of  the  radar  station,  set  NTRAP  = 
1.  The  measured  characteristics  of  the  layer  must  then  be 
entered  where, 

HI  =  The  base  height  of  the  layer  above  mean  sea  level 
(kilometers). 

SNI  =  The  value  of  refractivity  at  the  base  of  the  layer. 

DH  =  The  thickness  of  the  layer  (kilometers). 

DN  =  The  change  in  refractivity  through  the  layer  (always  a 
positive  quantity  for  the  inversion  layer). 


And  refractivity,  N,  is  defined  by 

„  77.  6  r  ,  4810  e  n 
N  =  ~  [p  +  — ^ ] 


where 


T  =  The  air  temperature  (deg.  Kelvin), 
p  =  The  air  pressure  (millibars), 
e  =  The  water  vapor  pressure  (millibars). 


90 


AFWTR  has  programs  already  available  to  calculate  N. 

The  local  earth  radius  used  is  seen  to  be  6378.  1 6 3  km  (page 
95  )  but  can  be  redefined  to  suit  local  geographical  latitudes. 

If  the  characteristics  of  the  inversion  layer  are  included  the 
program  will  indicate  if  the  radar  signal  can  be  trapped  and, 
if  so,  it  will  calculate  the  elevation  angle  at  and  below  which 
trapping  could  occur. 

If  no  layer  is  present,  or  if  it  is  desired  not  to  include  the 
effect  of  the  layer  set,  NTRAP  =  0. 

b.  The  second  option  assumes  that  there  is  no  direct  measure¬ 

ment  of  SNS  or  any  radiosonde  profile  information  available. 
In  this  case  there  must  be  a  value  of  refractivity ,  SNL,  ob¬ 
tained  at  the  radiosonde  launch  site.  This  value,  SNL,  is  set 
in  and  also  the  height  above  mean  sea  level,  HL  (km),  where 
SNL  is  measured.  Also,  of  course,  set,  NTRAP  =  0,  since 
no  inversion  layer  information  is  used. 

2 .  2  Pages  97  -  102,  Generalized  Mainline  Program  for  Polynomial 

Curve  Fit 


This  program  is  used  to  generate  the  values  of  four  constants  used  in 
the  real-time  running  program  (page  110).  These  constants,  CKO,  CK1,  CK2, 
and  CK3  are  functions  of  the  equivalent  radar  station  value  of  refractivity, 

SNS.  They  are  used  to  provide  an  optimum  calculation  of  refraction  errors 
for  short  tracking  ranges  and  small  elevation  angles.  These  constants  can 
be  obtained  from  graphs  or  tables  in  place  of  using  this  generating  program. 
The  program  is  included  herein  since  it  facilitates  direct  and  accurate  de¬ 
termination  of  these  constants  without  manual  intervention. 

2 . 3  Pages  103  -  108:  Complete  Magnetic  Tape  Routines 

This  program  is  designed  to  handle  all  magnetic  tape  operations.  It 
has  been  written  in  metasymbol  making  use  of  the  interlace  feature  and  the 
extended  mode  capability  but  does  not  use  or  affect  the  interrupt  system  in 
any  way.  It  contains  all  of  the  present  Fortran  capabilities  plus  several 
more  and  will  execute  all  magnetic  tape  operations  much  more  rapidly  and 
efficiently  than  do  the  Fortran  routines. 

It  is  used  here  to  output  the  necessary  parameters  calculated  in  the 
initialization  program. 
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THE  REFRACTION  CORRECTION  PROGRAM 


3.0 

3.  1  Pages  109  -  112;  SURC  Refraction  Correction  Program 

This  is  the  real-time  refraction  correction  program.  It  will  use 
constants  derived  from  the  Initialization  Program  (Section  2.  1)  and  the 
Polynomial  Program  (Section  2.2).  Given  the  uncorrected  radar  range,  R 
(feet)  and  the  uncorrected  elevation  angle,  (degrees),  it  will  generate  the 
refraction  induced  errors  and  provide  corrected  range  and  elevation  angle 
data.  It  will  also  generate  the  doppler  velocity  error  angle,  6  (DELTA) 
which  can  be  then  used  to  correct  doppler  or  range  rate  velocity  measure¬ 
ment.  (Reference:  Rowlandson,  L.  G,  ,  and  J.  R.  Herlihy,  "Refraction- 
Induced  Tracking  Errors  and  Correction  Methods  for  the  Air  Force  Western 
Test  Range,  "  December  1968,  ESD,  TR-69-52.) 

3 .  2  Pages  113  -  114  Converter  Program 

This  program  written  in  metasymbol  for  the  SDS-930  is  used  to 
convert  the  radar  data  supplied  on  AFWTR  tapes  (FPS  -  16  data,  36  bit 
floating  point  words)  to  a  format  compatible  with  the  input  requirements 
for  the  SDS-930  (48  bit  floating  point  words). 

3  3  Pages  115  -  120  Complete  Magnetic  Tape  Routines 

This  program  is  identical  with  2.  3.  It  is  used  to  input  from  magnetic 
tape  the  parameters  from  the  initialization  program  and  the  raw  radar  data. 
Under  sense  switch  control,  this  program  will  output  the  results  for  possible 
future  display  purposes. 
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COMMENTS 


4.  0 


The  above  program  is  being  used  in  the  SDS-930  to  analyze  AFWTR 
radar  tapes  at  a  data  rate  of  8  radar  data  points  per  second  (approximately 
124  ms  execution  time  per  data  point).  If  the  program  were  run  in  machine 
language  on  the  930  it  could  readily  handle  a  data  rate  of  10  points  per  second. 

We  believe  that  in  machine  language  on  the  IBM  computer  at  AFWTR 
the  program  should  execute  in  less  than  50  ms  because  of  the  reduced  times 
required  for  division  and  multiplication. 

This  SURC  refraction  correction  program  does  not  smooth  the  radar 
data  but  gives  a  direct  correction  for  refraction  on  a  one-to-one  relationship 
with  the  input  data. 

It  should  be  pointed  out  that  the  program  was  developed  for  use  in  the 
metric  system  requiring  that  conversions  be  made  from  feet  and  degrees  to 
kilometers  and  radians.  This  conversion  is  presently  carried  out  within  the 
program  thereby  making  it  directly  compatible  with  AFWTR  radar  data  for 
raw  radar  range  input  in  feet  and  raw  tracking  angle  in  degrees.  The  metric 
system  was  used  within  the  calculations  because  the  CRPL  exponential  profile 
data  and  ray-traced  error  data  were  metric. 

It  should  be  repeated  that  all  height  information  put  into  the  initializa¬ 
tion  program  (2.  1)  must  be  in  kilometers .  Of  course,  AFWTR  may  wish  to 
make  the  conversion  necessary  to  use  height  data  in  feet. 
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aUOB. 

^rewind  bo. 

FORTRAN  BO/LB. 

refraction  CORRECTION 
INITIALIZATION 
Q000l*01  »•>-  PEISEL 


THE  INITIALIZATION  PROGRAM  READS  9  DATA  CARDS 

plus  a  set  of  data  to  be  used  by  the  polynomial 

curve  fit  SUBROUTINE.  FROM  THESE  IT  CALCULATES 
SOME  NECESSARY  PARAMETERS  FOR  THE  ReFCOR 
PROGRAM  AND  OUTPUTS  THEM  ON  TAPE  UNIT  1 .  THE  9 
INPUT  PARAMETERS  TO  THJS  PROGRAM  ARE  IN  Fio.0 
FORMAT.  THE  UNITS  ARE  KILOMETERS  »• 

HS  •  STATION  HEIGHT  ABOVE  MEAN  SEA  LEVEL 
HL  *  RADIOSONDE  LAUNCH 

STATION  HEIGHT  ABOVE  SEA  LEVFL 
SNS  *  STATION  VALUE  OF  ReFRACTIVITY 

SNL  *  RADIOSONDE  LAUNCH 

STATION  VALUE  OF  REFRACTIVITY 

ntrap  •  l  IF  layer  IS 

TO  BE  INCLUDED/  0  IF  NOT 

hi  *  layer  height  above  mean  sea  level 
SN1  B  REFRACTVITY  VALUE  at  height  hi 
DH  -  THICKNESS  OF  THE  LAYER 

DN  a  CHANGE  IN  REFRACTIVITY 

through  the  layer 

C 

C 

DIMENSION  DATA(bO)/KXLG( lO)i YX(0/3) 

EQUIVALENCE  (DATA! i )/XXLG( 1 > >/ 

(DAT A ( J 1 ) / AnGL  > / <DATA(12)/SnS) 

EQUIVALENCE  ( DATA ( 1 3 ) / HS ) / ( DATA ( 14 ) / Hi ) / ( DATA ( lb ) / SNl ) 
EQUIVALENCE  ( DAT A ( 1 6 ) / DH ) / ( DATA ( 1 7 ) * DN > / (DATA ( jg ) / CC ) 
EQUIVALENCE  (DATA( 19) /CKO)/ (DATA (20)/ 

CKJ  )/ ( DA  T  A ( 2 1 ) / CK2 ) 

EQUIVALENCE  ( DAT A ( 22 )  /  CK3 )  /  (DATA (23)/ 

PI )/ (DATA(24)/PIPI ) 

EQUIVALENCE  ( D AT  A ( 25 ) / P I  2 ) / (PATA(26)/ 

DTOR ) / ( DAT  A ( 27 ) / RTOD ) 

EQUIVALENCE  ( D AT A ( 28 ) / ERAD I  US ) / 

(DATA (29) /SGI )/ ( DAT  A ( 30 ) / SQ2 ) 

EQUIVALENCE  ( DAT  A ( 3 1 ) / AE ) / (QATA<32)/ 

DH2)/ ( DATA ( 33 ) / AEHS ) 

EQUIVALENCE  (DATA ( 34 ) / AEHl ) / (DAT  A (35)/ 

FTOKM), ( DAT  A ( 36 ) / AKTOE ) 

EQUI VALLNCF  (DATA { 37 ) / EL  TRAP ) / 

( D  A  T  A ( 3  8 ) / N ) / (DAT A (39) /NTRAp) 


9-1 


J 
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c 


c 


c 


c 


PRINT  60 

60  F9RMAT(lHi,/,20X,$R&FRACTieN  C8RRECTI8N$, 
//23X,$INITIALIZaTIQN$, 

1///) 

READ  10/HS^HLiSNSiSNL>NTRAP/Hl/SNl/DH/DN 
10  F8RmAT(F10.0) 


P I  *3 • 14159265 
P I P I »2»  *P I 
Pl2»PI/2» 

DTSRbP 1/180# 

RT8D«180»/PI 
AKT0F«328O»8236 
FT8«M» 1  * /AKT0F 
ERADIUS.6378. 163 
DH2«DH+DH 
ANGl,«  1  • 

CCCfOt 

DNUOt 

ELTRAPbO. 

N*  lQ 
XN-N 

XXN-AU8GF(XN) 

D8  20  FI*2/N 

20  XXLG(FI )*1.-AU8GF(XN-FI+1. )/XXN 


I F ( SNL ) 32/ 31 / 32 

32  DNU-7.32*EXPF(  .005577*SNL) 
CCC-AL8GF(SNL/(SNL+DNn ) 
SNS*SNL*EXPF(-CCC#(HS-HU) ) 

31  DNNb"7.32#EXPF( .005577*SNS) 
CC»AL8GF(SNS/(SNS4.DNN)  ) 

SQ1«SQRTF(CC*ERADIUS*.5> 

SQ2bSQRTF(CC/(2.*ERADJUS) ) 

XK«1 •/ ( 1 •+ (ERADIUS/1 *0E  06 >M (SN1-SNS)/(H1-HS) ) ) 

AE»XK*ERADIUS 

AEHS-AE+HS 

AEH1.AE+H1 


IF(NTRAP)33/33/37 

37  IF(DN*l»0E-06-DH/ERADIUS)33/34#34 

33  PRINT  35 

35  FORMAT (5X/$N8  TRAP® ) 

G8  T8  36 

34  Ell»SQRTF( (2.*DN*1.0E-06)-(2.»DH/eRADIUS) ) 
EUTRAP«C8SF(ELl)/( 1#+(SNS.SNX >* 

1.0E*06-(H1-HS)/ERADIUS> 
EUTRAPbATANF (SQRTF ( l»^ELTRAp*ELTRAp)/ELTRAp) 


95 


PAGE  3 


PRINT  38/ELTRAP 

38  F0RMAT(5X,$TRAP$,5X,$ElEVATieN  „  $/Fl2.6) 

C 

36  CALL  P0LYFIT(SNS/VX) 

CK0»YX(0)«1.0E-02 

CK1pYX(1) 

CK2*YX  ( 2 )  *1  i0E>02 
CK3hYX(3)*1*0E-02 
c 

PRINT  61>SNS/SNL>HS/HL>AE/CC>DNN* 
CCC/0NI/CK0/CK1/CK2/CK3 
61  FORMAT (///5X>$SNS  *  $>Fl2»5,/> 

5Xi$SNL  «  $>F12«5//^5X/$HS  «  $/ 

1F12.5///5X/$HL  *  *^F12.5///5X/ 

$AE  *  $*F12.5*/j5X,$CC  ■  $, 

2F12»5//>5X/$DNN  ■  $/ 

3F12.5>/>5X^$CCC  ■  $>F12t5//>5X* 

$DNI  -  $/Fl2»5>/.»5Xi$CK0  *  $* 
4F12.5>//5X>$CK1  ■  $iFl2.5///5X/ 

$CK2  «  $/F12t5//^5X^$CK3  «  $,F12.5) 
C 

CALL  WBT(DATA(1),100>2) 

C 

CALL  ENDF ( 2 ) 

CALL  REW ( 2 ) 

CALL  EXIT 
END 
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SUBRBUT I NE  POLYF I T ( SNS;  YX ) 

GENERALIZED  MAINLINE  PROGRAM 
FOR  'POLYNOMIAL  CURVE  FIT' 

DIMENSION  YX ( 0/3 ) 

DIMENSION  X(10Oj;YUO0);A(ll>, 

KTITlE(10);XX(20);YY(20); ITAB<3) 

LL«0 


FORMAT  STATEMENTS 


1  F  BRMAT ( 1QA4 ) 

2  F BRMAT ( 2F 1 0 • 0 ) 

3  FORMAT (4F1Q.0) 

4  FORMAT ( //// ( $  A  (  $;  1 2;  * )  •  $;E15*8  )) 

5  F0RmaT(1H1;4QX;10A4///V$  THE 

COMPLETE  SET  OF  DATA  INCLUDES*; 15;* 

1  POINTS.*; 10X;*X$;24X;$Y$//(45X;F20.8;5X;F20*8) ) 

6  FORMAT ( //V*  XM I N  ■$; F20 . 8; 20X; 

$XMAX  •$;F20t8/$  YMIN  «*;F20*8;20X 


1;*YMAX 
12  FBRMAT( 


*; F20 • 8/ / / ) 

1H1$  DATA  P0INT  NUMBER*; 

1 4; *  AND  THE  NEXT*; J5;$ 

ITIVe  points  will  BE  used  to  find  the 
14  FORMAT ( //$  XTEST  «  $;El5.8/$ 

YTEST  ■  $; E 1 5 • 8/$  CALX 
16  FORMAT (///////*  THE  FOLLOWING 

ARE  THE  COEFFICIENTS  OF  A  POLVNOMIA 
1L  OF  DEGREE*; 13;*  FOUND  AS  A 


CONSECU 
CURVE  FIT.*) 

*  $; E 15 • 8  ) 


CURVE  FIT  TO  THE  DATA  POINTS  GIVEN. $ 


2) 

19  FORMAT (1H1;*#*#**#  THE  CODE  SPECIFIED 
BY  THE  USER  IS  NOT  PERMISSJ 

1BLE.$//1H1 ) 

26  FORMAT ( ////*  THESE  POINTS  ARE*; 

10X;$X*;24X;$Y$;//(15X;F20.8;5X; 

1F20.8) ) 

55  FORMAT ( 1615 ) 

64  FORMAT ( /////*  THE  FOLLOWING  CASES 

PROVIDE  A  CHECK  for  this  pblynb 
1MIAL  CURVE  FIT.*) 

67  FORMAT ( ///$  **********  X  TEST 

VALUE  eQUALING$;e15.8;$  IS  BUTSID 
IE  THE  RANGE  OF  INpUT  DATA;  THUS 

NORMALIZING  IS  I NC0RRECT . $/$  THE  R 
2ANGE  OF  THE  X  VALUES  OF  THE  INPUT 

DATA  I S*; E 1 5 . 8; $  T0$;E15.8) 

I T  Ag ( 1 ) »4HSAME 
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I  T A 0 ( 2 ) ?4HNE W 
I  TAB ( 3 ) *4HEND 
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READ  TOTAL  NUMBER  OF  POINTS  IN  COMPLETE  DATA  SET. 

read  problem  title, 
read  complete  set  of  data 

POINTS  ( X/ Y )  to  be  used, 
data  POINTS  1  THROUGH  t  NUMBER  *  MUST  BE  READ  IN 
EVEN  IF  'L*/  THE  FIRST 

DATA  POINT  TO  BE  USED/  IS  GREATER  THAN  1. 

EACH  DATA  CARD  CONTAINS  ONE 

POINT  THUS  HAS  AN  X  AND  A  Y, 

DATA  SHOULD  BE  ARRANGED  IN  ASCENDING  ORDER. 

THE  CARD  FORMAT  IS  {2F1Q.Q). 

9  READ  1/KTITLE 
READ  55/NUMBER 

READ  2/ (X( I ) ,  Y  ( I )/ I  *1, NUMBER) 

PRINT  5/KTITLE/NUMBER, ( X ( I ) , Y ( I ) , I  *  1 # NUMBER ) 


read  highest  degree  of  fit 

DESIRED/  THE  NUMBER  OF  DATA  POINTS 
TO  Be  used  IN  CURVE  FITTING, 
the  SUBSCRIPT  OF  the  FIRST 
data  point  to  be  used,  and  the 
number  of  test  cases, 
read  xmin,  xmax,  ymjn,  ymax  of  data  set. 

20  READ  55,MM,L,N,NCHECKS 

READ  3/XMIN, XMAX, YMJN, YMAX 
NN"N- 1 

PRINT  12/L/NN 

PRINT  26/  (X(I  ),YU  )  /  I  ■L/L  +  NN  ) 

PRINT  6/XMIN/XMAX/ YMIN/YMAX 


NORMALIZE  X  AND  Y  VALUES  OF  DATA  POINTS 

DO  101  I»L/L+NN 
X (  I  )  « ( X ( I ) -XMIN)/(XMAX"XMJN) 

101  Y< I )= ( Y< I ) "YMlN)/( YMAX*YMIN) 
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C  READ  IN  ALL  TEST  P0 1  NTS  WHICH 

WILL  BE  USED  T0  EVALUATE  THE 
C  RESULTING  POLYNOMIAL* 

C 

IF(NCHECKS)60*60*6l 
61  XX(1)pSNS 
YY(1)«0.0 


fit  a  polynomial  of  degree 

M  THROUGH  the  set  of  data  pbints 

60  DO  21  M»2* MM 
DO  62  1*1*11 
62  A ( I } *0 • 

PRINT  16*  M 

CALL  LSQ(M,N,L*X,Y,A) 

PRINT  4* (I  *  A ( I) *  I  *  1  *  M+l ) 

IF (NCHECKS)21*21* 11 
11  PRINT  64 

11  continue 

DO  100  I  s  1  * NCHECKS 


NORMALIZE  TEST  VALUE  8F  XX 

TO  FIT  WITHIN  DATA  RANGE  AND 
C  CALCULATE  F(XX) 

C 

CALX»  0 • 

XTEST« (XX( I )-XMIN)/(XMAX-XMlN) 

IF (XTEST)65>66>66 
66  IF (XTEST-i • )68*68*65 
65  TYPE  67*XX( I ),XMIN,XMAX 
GO  TO  100 
68  DO  63  J* 1 *  M+ 1 
63  CALX»A(J)*XTEST*MJ"1  )+CALX 
CALX«CALX*(YMAX-YMIN)+YMIN 
C  PRINT  14*  XX ( I ) *  YY ( I ) *  CAL* 
IF(M-3)401*400*401 
C  400  TYPE  401*LL#CALX 

400  CONTINUE 
YX(LL)«CALX 

401  F8RMAT(2X>$K$, I2>$  *  $,El5*8) 

100  continue 
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21  continue 

LL=LL+1 


READ  CODE  CARD,  THE  OPTIONS  INCLUDE; 
SAME  DATA  SET 
NEW  DATA  SET 
END 

READ  It  I  CODE 
DO  17  1*1/3 

I F ( I  CODE- I  TAB ( I ) ) 17/ 18/ 17 

17  CONTINUE 
PRINT  19 
GO  TO  1113 

18  IF( 1-2)20/9/ 1113 
1113  RETURN 

END 
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SUBROUTINE  LSQ(MiN/LiXiY/A) 

LEAST  SQUARES  POLYNOMIAL  CURVE  FIT 

FITS  AN  m-th  degree  polynomial 
through  N  points 
starting  WITH  the  l-th 

dimension  X(l)iY(l)iA(l)iBUl*18) 

clear  the  matrix 


DO  2  InliM+1 
DO  2  J-l/M+2 

2  B( It J)»0 

FORM  THE  LOWER  MATRIX 

DO  4  K«LiL+N-l 
XI -1 

DO  4  I • 1 j  M  +  l 
XJ*1 

DO  3  J»1>I 
BUiJ)«B(  U  J)+Xl*XJ 

3  X J*X J*X ( K ) 

B{ IjM+2)»B( IiM+2)+Xl*Y(K) 

4  X  I  «X I *X ( K ) 

MAKE  MATRIX  SYMMETRIC  AND  SOLVE 

DO  5  I  ■ 1 j  M+ 1 
DO  5  J» 1/ I »1 

5  B ( Ji I ) *B ( I / J ) 

CALL  S0LVELSQ(M+liBi INDlC) 

I F ( INDIC)7*7j8 

MATRIX  IS  SINGULAR 

8  TYPE  50 

50  FORMAT (//$N0  S0LN,  EQS  SINGULAR*//) 
RETURN 

7  DO  6  I » li M+l 

6  A ( I ) »B ( I  4  M  +  2 ) 

RETURN 

END 
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SUBROUTINE  S6LVEISQ(NA, A/ INDIC) 

modified  version  of  solve 

FOR  USE  WITH  LSQ 

DIMENSION  A ( 1 1/ 12 ) 

DO  205  I*1/NA 
I F ( A ( I / I ) )200/2Q6/2Q0 

200  Pp 1/A (  1/  I  ) 

D8  201  Us  I  / NA+ 1 

201  A( I, J)»P*A( 1/ J) 

DO  202  K-l/NA 

I F ( I*K)2Q3/202/203 

203  Q* A ( K/  I  ) 

DO  204  J* I / NA+1 

204  A(K/j)»A(K/J)«Q*A( 1/ J) 

202  CONTINUE 

205  CONTINUE 
I ND I C  »0 
RETURN 

206  IF(NA-I )3Q0/304/3Q0 

300  DO  303  Ls 1+1/ NA 

I F ( A ( L/ I ) ) 30 1 / 303/ 301 

301  DO  302  J« I / NA+ 1 
HOLD" A ( I / J ) 

A  (  I  ,  J  )  «  A  (  L/  U  ) 

302  A(L,J)“H0LD 
GO  TO  200 

303  CONTINUE 

304  I ND I C* 1 
RETURN 
END 
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aEOF. 

AMETA920  SI/LO/BO. 

*  C8MPLETE  MAGNETIC  TAPE  ROUTINES 

*  PEISEL 

« 

*  THESE  ROUTINES  ARE  capable  0P 

transmitting  Data  td 

*  AND  from  MAGNETIC  TAPE  IN  EITHER 

A  BINARY  OR  DECIMAL  MODE. 

*  any  number  of  records  or  files 

MAY  BE  SKIPPED  in  EITHER 

*  DIRECTION/  THE  TAPE  UNIT  MAY  Be 

REWOUND/  AN  END-OF-FILE 

*  may  be  detected  or  written  IF  desired. 

IN  ADDITION/  IF 

*  THE  WORD-COUNT  IN  the  RECORD  IS 

UNKNOWN/  THE  ENTIRE  RECORD 

*  may  Be  read  and  the  number  of  WORDS 

RETURNED  IN  K8UNT 

*  BY  ENTERING  the  ROUTINE  WITH  kount 

■  -1.  AN  INDICATOR 

*  Is  iNCeRPoRATED  Te  ALLeW  REC0VERY 

from  TRANSMISSION  ERRORS. 

*  THE  FOLLOWING  IS  A  LIST  OF  TH£ 

call  statements  as 

*  THEY  SHOULD  APPEAR  IN  THE  FORTRAN 

program  with  the 

*  ARGUMENTS  DEFINED: 

*  data  ■  the  dimensioned  data  array 

*  KOUNT  «  WeRD  count 

*  MaGt  ■  LOGICAL  MAGNETIC  TAPE  unit 

*  KFILE  ■  0  NORMAL  CORRECT  TRANSMISSION 

*  «  -1  READ  ERRSR 

*  «  +1  END-OF-FILE 

«  N  .  +  OR  -  THE  NO.  OF  RECORDS 

OR  FILES  TO  BE  SKIPPED 

*  CALL  RBT (DATA(I)/KOUNT/MAGT/ KFILE) 

*  CALL  RDT(DATA( 1 )/KOUNT/MAGT/KFlLE) 

*  CALL  WBT(DATA( 1 )/K0UNT/MAGT) 

*  call  wdt<dataU)/Kount/Magt> 

*  CALL  SKF(MAGT*N) 

*  CALL  SKRC(MAGTzN) 

*  CALL  ENDF(MAGT) 

*  CALL  REW(MAGT) 


RORG  0 

TFT  OPD  QA013610 
$WBT  PZE 

skr  MEM 


WRITE  TAPE  BINARY  ENTRY 
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$WDT 

PZE 

WRITE  TAPE  DECIMAL  ENTRY 

skr 

mem 

$RBT 

PZe 

Read  tape  binary  entry 

skr 

mem 

$RDT 

pze 

Read  tape  decimal  entry 

LDa 

erct 

maximum 

number  of  transmission  errors 

STA 

ERRCt 

ERROR  COUNT 

UDa 

Mem 

ADD 

MeM 

MULTIPLY  BY  2 

RCH 

0402 

0  INTO  B/  a  into  X 

STB 

MEM 

restore  mem  to  zero 

ADD 

three 

read  or  write  switch 

STA 

cntrl 

set  switch 

LDA 

T0P12/2 

retrieve 

PR0PeR  TAPE  OPERATION 

STA 

rwt 

LDa 

RDT/2 

obtain  return  address 

STA 

RBT 

LDX 

E8ADR 

address  -  erasable  core 

SKN 

CNTRL 

STB 

*3/2 

clear  KFILE  INDICATOR 

LDa 

*1/2 

2ND  ARGUMENT 

STA 

RC0UNT 

W0RD  COUNT 

ETR 

MAXWC 

CLEAR  NEGATIVE  WORD  COUNT 

RSH 

10 

EXTENDED  word  count 

mrg 

TQPll 

EXTENDED  word  count-eom 

STA 

M0DE 

STB 

PW0RD 

LDa 

*E0 ADR 

ADDRESS  »  1ST  ARGUMENT 

ETR 

ADRMSK 

MRG 

PW0RD 

NORMAL  W0RD  COUNT 

STA 

PW0RD 

INTERPLACE  POT  word 

LDA 

*2/2 

3RD  ARGUMENT 

STA 

magt 

LOGICAL  MAGTAPE  UNIT 

MRG 

RWT 

merge  with  tape  operation 

STA 

RWT 

SKN 

CNTRL 

read  -  write  switch 

BRu 

rerun 

LDa 

T0P1O 

beginning  of  tape  test 

mrg 

magt 

STA 

$  +  2 

BRti 

delay 

N0P 

TAPE  OPERATION 

BRM 

eras 

BRU 

RERUN 

ERAS 

P^E 

ERASE  TAPE  forward 

LDA 

T0P9 

MRG 

magt 
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STA 

$+1 

N8P 

TAPE  OPERATION 

POT 

dummy 

brr 

ERAS 

RERUN 

brm 

DELAY 

rwt 

NOP 

TAPE  OPERATION 

M0DE 

N0P 

EXTENDED  MODE  EDM 

P0T 

PW0RD 

lda 

T0P7 

MRG 

MAGT 

STA 

$  +  1 

N8P 

TAPE  READY  TEST 

BRU 

$+3 

CZT 

0 

CHANNEL  ZERO  TEST 

BRU 

$  3 

TFT 

0 

TAPE  END-OF-FILE  TEST 

BRu 

file 

CET 

CHANNEL  ERROR  TEST 

BRU 

TRERR 

SKN 

CNTRL 

READ  •  WRITE  SWITCH 

SKN 

RC0UNT 

BRR 

RBT 

ASC 

PIN 

RC0UNT 

RETRIEVE  LAST  ADDRESS 

LDA 

*E0ADR 

ETR 

adrmsk 

SUB 

RC0UNT 

obtain  actual  word  count 

cna 

STA 

*1/3 

RETURN 

COUNT 

TO  MAIN  PROGRAM 

BRR 

rbt 

RETURN 

file 

STA 

*3/3 

SET  END-OF-FILE  INDICATOR 

BRR 

RBT 

RETURN 

trerr 

brm 

DELAY 

transmission 

error 

ROUTINE 

LDa 

MAGT 

SKN 

CNTRL 

READ  -  WRITE  SWITCH 

BRU 

SCN 

MRG 

T8P3 

erase  reverse  tape 

STA 

$+1 

N0P 

TAPE  OPERATION 

EXU 

MODE 

extended  mode  eom 

P0T 

PW0RD 

brm 

ERCNT 

BRU 

rerun 

RETURN  TO  REWRITE 

SCN 

MRG 

T0P4 

SCAN  REVERSE 

STA 

$  +  1 

N0P 

TAPE  OPERATION 

P0T 

ERR 
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BRM 

ERCNT 

BRU 

rerun 

ERCNT 

PZE 

min 

ERRCT 

skn 

errct 

BRU 

$4-2 

BRR 

ERCNT 

SKN 

CNTRL 

BRU 

ERRD 

lda 

magt 

MRG 

T0P9 

STA 

$  +  2 

brm 

NOP 

DELAY 

EXU 

MODE 

POT 

PWORD 

lda 

ERCT 

STA 

errct 

BRR 

ERCNT 

ERRD 

lda 

ONES 

BRU 

FILE 

$SKT 

RES 

0 

$SKF 

PZE 

lda 

ONES 

LDB 

$-2 

STB 

$  +  2 

BRU 

$  +  3 

$SKRC 

PZE 

CLA 

STA 

CNTRL 

LDX 

E0ADR 

lda 

*0/  2 

STA 

MAGT 

MRG 

T0P6 

STA 

SCAN 

STA 

RDF 

LDa 

*1/2 

STA 

rcount 

skn 

RCQUNT 

BRU 

CAX 

$  +  6 

lda 

T0P4 

MRG 

magt 

STA 

SCAN 

BRU 

$  +  5 

skg 

ZERO 

BRR 

CNA 

CAX 

SKRC 

RETURN  TO  REREAD 
ERROR  COUNT  R8UTINE 

ARE  ERROR  TRIALS  DONE 

TRY  AGAIN 

failed  after  10  tries 


ERASE  TAPE  F0RWARD 


ERASE  passed  BAD  SP0T 


restore  error  count 

CONTINUE 

RETURN  WITH  KFJLE  =  -1 


file  switch 

obtain  return  address 


RECORD  SWITCH 
SET  SWITCH 

ADDRESS  -  ERASABLE  core 
1ST  ARGUMENT 
LOGICAL  mag  TAPE  UNIT 
SCAN  FORWARD 


2ND  ARGUMENT 

record  or  file  count 
word  Count 


SCAN  REVERSE 


ZERO  COUNT  RETURN 
SKIP  forward 


106 


PAGE  5 


BRM 

DELAY 

SCAN 

NOP 

TAPE  OPERATION 

P0T 

ERR 

skn 

CNTRL 

record  0R  file  switch 

BRU 

SKR 

brm 

DELAY 

TFT 

0 

TAPE  END-OF-FILE  TEST 

BRU 

$+2 

BRU 

SCAN 

BRX 

SCAN 

skn 

RC0UNT 

RECORD  OR  FILE  SWITCH 

BRR 

SKRC 

RETURN 

RDF 

N0P 

PLACE  READ- 

head  beyond  file  mark 

P0T 

ERR 

BRR 

SKRC 

RETURN 

SKR 

BRX 

SCAN« 1 

BRR 

SKRC 

RETURN 

$rew 

PZE 

rewind  entry 

lda 

*E0  I ND 

STA 

magt 

MRG 

T0P5 

REWIND 

STA 

$  +  2 

BRM 

DELAY 

N0P 

REWIND  TAPE 

BRR 

REW 

RETURN 

$ENDF 

PZE 

WRITE  END-OF-FILE  ENTRY 

LDA 

*EOIND 

STA 

magt 

MRG 

T0P1O 

beginning  of  tape  test 

STA 

$  +  1 

N0P 

TAPE  OPERATION 

BRM 

ERAS 

lda 

magt 

MRG 

T0P8 

write  tape  decimal 

STA 

$+2 

BRM 

DELAY 

N0P 

write  END-OF-FILE 

P0T 

TMC 

BRR 

ENDF 

RETURN 

delay 

PZE 

TEST  STATUS  OF  TAPE 

LDA 

T0P7 

mrg 

MAGT 

STA 

$<♦•1 

N0P 

TAPE  READY  TEST 

CAT 

CHANNEL  ACTIVE  TEST 

BRU 

$-2 

BRR 

DELAY 

INTER 

FORM 

10/ 14 

interlace  format 
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PWQRD 

INTER 

0/0 

P0T  WORD 

ERR 

inter 

1 ,  err  +  i 

SCAN  P0T  W0RD 

RES 

1 

INPUT  1ST  W0RD  FROM  SCAN 

DUMMY 

inter 

600/0 

ERASE  BEY0ND  LOAD  POINT 

TMC 

INTER 

1/TMOl 

data 

017000000 

E0ADR 

eQu 

071 

E0IND 

EQU 

074 

adrmsk 

EQU 

027 

ADDRESS  MASK 

ZERO 

EQU 

023 

ONES 

EQU 

026 

THREE 

EQU 

0342 

magt 

RES 

1 

cntrl 

RES 

1 

RC0UNT 

RES 

1 

ERRCT 

RES 

1 

MEM 

data 

0 

ERCT 

data 

-10 

MAXWC 

data 

077777 

T0P1 

WT8 

*0/0/4 

T0P3 

ERT 

*0/0/4 

T0P13 

WTD 

*0/0,4 

T0P4 

SRB 

*0/ 0, 4 

T0P2 

rtb 

*0/0/4 

T0P5 

REW 

Q*0 

T0P12 

RTD 

*0/0/4 

general  tape  operations 

T0P6 

sfb 

*0/0/4 

T0P7 

TRT 

0/0 

T0P8 

WTD 

*0/0/ 1 

T0P9 

eft 

*0/0/4 

T0P1O 

btt 

0/0 

T0P11 

EOM 

014000 

m  m  f»  * 

end 

aEOF. 

^rewind  bq. 

aF0RTL6AD  MAP, IMAP, BI. 
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A1J8B  • 

^REWIND  B0. 

FORTRAN  80/L0. 

refraction  CORRECTION 
REFCOR 

QOOOl-Ol  --  PEISEL 


the  refcor  program  Reads  the  parameter  list 

OUTPUT  FROM  the  INITIALIZATI9N  PROGRAM  FROM  TAPE 
UNIT  1  AND  the  MEASURED  ELEVATI8N  AND  RANGE  DATA 
FROM  TAPE  UNIT  2.  IT  TH£N  CALCULATES  THE  ERR8R 
AND  PRINTS  THE  RESULTS  IF  SENSE  SWITCH  4  IS  SET 
BR  OUTPUTS  THE  RESULTS  ON  TAPE  UNIT  3  IF  SENSE 
SWITCH  3  IS  SET. 


DIMENSION  DATA(5Q)/XXLG( 10),LDATA( 12) /KDATA (34) 
EQUIVALENCE  (DATA(1)/XXLG(1) >/ 

(DATA(U)/ANGL)MDATA(12)/SnS) 

EQUIVALENCE  (DATA(13)>HS), (DATA(14)#Hl), (DATA(15)/SN1) 
EQUIVALENCE  ( DATA ( 16 ) / DH ) / ( DATA ( 17 ) / DN ) / ( DATA ( 18 ) / CC ) 
EQUIVALENCE  ( DATA ( I9 )/ CKO ) / ( DATA ( 20 ) / 

CK1)/ <DATA(21)/CK2) 

EQUIVALENCE  ( DAT A ( 22 ) / CK3 ) /  (DATA (23)/ 

P I  )  /  ( DAT  A ( 24 ) *  P I P  I  ) 

EQUIVALENCE  ( DATA C25 ) / PI 2 ) /  ( DATA ( 26 ) / 

DTOR ) / (DATA(27)/RT0D) 

EQUIVALENCE  (DATA(28)/ERADIUS), 

( DATA ( 29 ) / SQ1 ) / ( DATA ( 30 ) *  SQ2 ) 

EQUIVALENCE  ( DATA ( 31 ) / AE ) /  (DATA (32)/ 

DH2)/ (DATA(33)/AEH$) 

EQUIVALENCE  (DATA(34)/AEH1 )/ (DATA (35)/ 

FTOKM)/ (DATA(36)/AKT0F) 

EQUIVALENCE  (DATA(37)/ELTRAP)/ 

(DATA(38)/N)/ (DATA(39)/NTRAP) 

EQUIVALENCE  ( KDATA ( 4 ) / T I  ME ) / ( KDAT  A ( 6 ) / 

VRANGE)/ (KDATA( 10J/VEL) 

EQUIVALENCE  ( KDATA ( 12 )/ VCR )/( KDATA ( 16 ) / 

VCEL)/ (KDATA( 18J/VDR) 

EQUIVALENCE  ( KDATA ( 20 ) / VDEL ) / ( KDAT A ( 22 ) / 

CRANGE)/ (KDATA(24)/CELEV) 

EQUIVALENCE  ( KDATA ( 26 )/ ORANGE ) / 

(KDATA(28 ) /DELEV)/ (KDATA (30) /RD IFF) 
EQUIVALENCE  ( KDATA ( 32 )/ ED  I FF )/( KDATA ( 34 )/ DELTA ) 

C 

PRINT  66 

66  FORMAT ( 1H1////48X# $CBMPaR I SON$/ 

//36X/$vandanBeRg  method  and  surc 
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1  METH6D$,/#45X,$F0R  CORRECT  1 ON 

0F$//i35X/$RANGE  AND  ELEVATION  REF 
SRaCj I  ON  ERRORS////) 

PRINT  68 

68  FORMAT (// 31 X/  $-•■*  VANDANBERG 

---$,8X, $-•-$, 6*/ $$URC$/7X, 
4$---$/ 8X#  $SURC  -  VAN$, 12X,$ANGLE$, 

6X#  $EL APSE$/ //  15X/$MEASURED$/ 
58X,$C0KRCCTED$i8X,$ERR0R$,8X,$CeRRECTED$, 
9X,$ERR0R$/8X# 

6$DIFFERENCE$#12X#$DELTA$#7X#$TIME$#///) 

LNS.-13 

CALL  RBT(DATA(l),lOQ/8/KFILE) 

IF ( KF ILE ) 10/ 12/ 30 

12  CALL  RBT(LDATA( 1 )# 12/ l^KFILE) 

I F { KF ILE) 10/20/30 

10  TYPE  11 

11  FORMAT ($READ  ERRORS) 

G0  TO  30 

20  CALL  C6NV(LDATA( 1 )/KDATA(  1 )  ) 

ELCV»VEL*DT0R 
RANGE«VRANGE*FTOKM 
S I  Ne,L*S  I NF  (  ELE  V  ) 

SG)3«EXPF(-CK1*ELEV) 

IF(eLEV-ANGL)53/53/54 

54  SINGAM=SINEL 
L0P  b  l 

GO  TO  5b 

53  SINGAM»SINEL+(CK0+CK2*EXPF(.CK3*RANGE) )*S03 
LQPbO 

55  COSGaN=SGRtP ( 1 . . ( S I NGAM*S I NGAM ) ) 
TANGAMbSINGAM/COSGAM 
gg=sqi*tangam 

TAUF»SNS*1 »0E-06/T  ANGAM 

SK*RaNGE*SQ2*C0SGaM 

gk»gg+sk 

DRANGE»SNS#1.0E-06/(CC*SINGAM)« 

(1*»GG*EXPF{-SK*(2.*GG*SK) )/GK) 

KORb-1 

IF (NTRAP J56/56/ 57 
57  ATP.ELEV 

IF(ELEV-ELTRAP)47/48/48 

47  ELEV-ELTRAP 

48  ALPHA *CBbF(ELEV)*AEHS/AEHl 

ALPHA ■ AT ANF ( ALPHA/SQRTF ( l.-ALPHA*ALPHA) ) 
el  i  =p 1 2" alpha 


no 


PAGE  3 


EL2.SQRTF(ELl*Ell»DN#2.QE-06+DH2/ERADIUS) 

RC»AEH1*C0SF( ALPHA )"AEHS#S INF (EUEV) 

TM*EU1/EL2 

ELEV«ATP 

K0P»O 

56  SINDEL^O. 

DROwO • 

D0  40  I  »2/  N 
DR*RANGE»XXLG( I ) 

DDR  *DR-DRO 
DRO»DR 

DR*DR-DDR/2 • 

IF(L0P>44/44/45 

44  SINgAM«SINEL-MCK0+CK2*EXPF(«CK3#DR)  )*SQ3 
CBSGAMbSQRTF( 1-SINGaM*SINGAM) 
TANGaM»SINGAM/C0SGAM 

gg*sqi*tangam 

TAUF*SNS*1»0E "06/TANG AM 

45  SK»DR*SQ2*C0SGAM 
GK«GG+SK 

TAUrTAUF*( 1»-GG*EXPF(-SK*(SK+2.#GG) )/GK) 
IF(K0P)4l/42/43 

42  lF(DR"RC)41/43/43 

43  SINDEU«SINDEL4SINF(ELEV-TAU*TM)*DDR 
K0P.1 

G0  T0  40 

41  SINDEL«SINDEL4-SINF(ELEV-TAU)*DDR 
40  CONTINUE 

VDR«VRANGE-VCR 

CRANGE*(RANGE«DRANGE)*AKT0F 

DRANGE«DRANGE*AKT0F 

RDIFF«VCR-CRANGE 

celev«sindel/range 

CELEV«ATANF(CEUEV/SQRTF( 1-CELEV#CELEV) ) 

deiev*elev-celev 

DELTA*  (TAU-DEt,EV)*RT0D 
DELEV*DElEV*RT0D 
CELEV*CELEV*RTOD 
VDEL-VEL-VCEL 
EDIFF«VCEL-CELEV 
I F ( SENSE  SWITCH  4)59/28 
59  I F ( LNS ) 58/ 27/  27 
27  PRINT  69 
69  F0RMAT(1H1) 

PRINT  68 
LNS, -16 

58  PRINT  61/VRANGE/VCR/VDR/CRANGE/ 

DRANGE/RDIFF,DELTA,TIME/VEL/VCEL/ 
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1VDEL>CLLLV> DELEV/ ED  IFF 
61  F0RMAT ( 2X,$RANGE$/Fi8t8/Fl7.8* 

F13.8>Fl7»8>Fl4*8>Fl5*8'F20.8> 

1F11.  W/*ELEVATI8N$/Fl6.8/Fl7t8/ 

Fl3»8/F17*8/F14»8/Fl5»8//) 

LNSsLNS+1 

28  IF(SENSE  SWITCH  3)31/12 

31  CALL  W6T(KDATA( 1 )/34/3) 

GO  TO  12 

30  CALL  REW(l) 

CALL  REW(2) 

IF (SENSE  SWITCH  3)32,33 

32  CALL  ENDF ( 3 ) 

CALL  RE W ( 3 ) 

33  CALL  EXIT 
END 
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aE8F. 

AMETA92Q  SI*LO*BO. 

*  converter 

*  72  BIT  FLOATING  POINT  WORDS  T8 

*  48  BIT  FLOATING  POINT  WORDS 

*  LDATA  TO  KDATA 


# 

RORG 

0 

$CBNV 

pze 

LDa 

eoadr 

ADDRESS  .  ERASEABLE  CORE 

STA 

KDATA 

LDa 

*KDAT  A 

STA 

LDATA 

ADDRESS  •  LDATA 

MIN 

KDATA 

LDa 

*KDAT  A 

STA 

KDATA 

ADDRESS  -  KDATA 

LDA 

#4 

STA 

WCNT 

4  TIMES 

LOOP 

SKR 

WCNT 

BRg 

$+2 

brr 

CONV 

RETURN 

LDX 

LDATA 

PHASE  1 

LDA 

0*2 

LRSH 

15 

ETR 

•  0577 

STA 

SAVE 

EXPONENT 

LDA 

1*2 

ETR 

=077770000 

CAB 

LDa 

0*2 

ETR 

■00077777 

lsh 

8 

mantissa 

LDX 

SAVE 

lde 

exponent  into  b 

LDX 

KDATA 

STB 

0*2 

STORE  FLOATING  POINT  WORD 

JL 

STA 

1*2 

in  reverse  order 

LDX 

LDATA 

PHASE  2 

LDA 

1*2 

LRSH 

3 

ETR 

■0577 

STA 

SAVE 

EXPONENT 

LDA 

1*2 

ldb 

2*2 

lsh 

20 

ETR 

■037777777 

MANTISSA 

LDX 

SAVE 

lde 

EXPONENT  INTO  B 
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LQX 

KDATA 

STB 

2*2 

store  floating  point  WORD 

sta 

3*2 

IN  REVERSE  order 

lda 

«3 

ADM 

LDATA 

MOVE  LDATA  POINTER 

LDA 

*4 

ADM 

KDATA 

MOVE  KDATA  POINTER 

# 

BRU 

LOOP 

# 

SAVE 

RES 

1 

LDATA 

RES 

1 

KDATA 

RES 

1 

WCNT 

RES 

X 

E0ADR 

Eug 

071 

END 
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aE0F. 

aMETA920  SI/L0/BB* 


* 

COMPLETE  MAGNETIC  tape  ROUTINES 

* 

# 

PEISEL 

# 

these  routines  are  CAPABLE  OF 

transmitting  data  to 

* 

AND  FROM  MAGNETIC  TAPE  IN  EITHER 

a  binary  or  decimal  mode. 

# 

ANY  NUMBER  OF  RECORDS  OR  FILES 

MAY  BE  SKIPPED  IN  EITHER 

# 

DIRECTION/  THE  TAPE  UNIT  MAY  Be 

REWOUND/  AN  END-0F-FILE 

* 

MAY  be  DETECTED  0R  WRITTEN  if  desired, 

IN  ADDITION/  IF 

* 

THE  W0RD.COUNT  IN  THE  RECORD  IS 

UNKNOWN/  THE  ENTIRE  RECORD 

# 

may  be  Read  and  the  number  of  words 

RETURNED  IN  K0UNT 

# 

BY  ENTERING  THE  ROUTINE  WITH  KOUNT 

«  -1.  AN  INDICATOR 

* 

IS  INCORPORATED  TO  ALLOW  RECOVERY 

FROM  TRANSMISSION  ERRORS, 

* 

THE  FOLLOWING  IS  A  LIST  OF  THE 

CALL  STATEMENTS  AS 

* 

THEY  SHOULD  APPEAR  IN  THE  FORTRAN 

PROGRAM  WITH  THE 

* 

ARGUMENTS  DEFINED; 

# 

data  *  the  dimensioned  data  array 

* 

KOUNT  »  WORD  COUNT 

# 

MAGt  «  LOGICAL  MAGNETIC  TAPE  unit 

# 

KFILE  «  0  NORMAL  CORRECT  TRANSMISSION 

# 

«  .1  READ  ERROR 

* 

■  +i  end-of "File 

# 

N  •  +  OR  -  THE  NO,  OF  RECORDS 

OR  FILES  TO  BE  SKIPPED 

# 

CALL  RBT(DATA( D/KOUNT/MAGT/KFILE) 

* 

call  rdt ( data (i)/Kount/Magt/ kfile) 

# 

call  wbtidatai i )/KounT/Magt) 

# 

call  wdtidatai i J/Kount/Magt) 

* 

CALL  SKF(MAGTzN) 

# 

Call  SKRC(MAGTzN) 

* 

CALL  ENDF(MAGT) 

* 

call  REW(MAGT) 

R0RG 

'  0 

TFT 

0PD 

04013610 

$WBT 

PZE 

WRITE  TAPE  BINARY  ENTRY 

SKR 

mem 

1 15 
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$WDT 

PZE 

write  tape  decimal  entry 

SKR 

MEM 

$RBT 

PZE 

read  tape  binary  entry 

skr 

mem 

$RDT 

PZe 

read  tape  decimal  entry 

lda 

ERCT 

maximum 

NUMBER  0F  TRANSMISSI9N  ERRORS 

STA 

errct 

error  count 

LDA 

mem 

add 

Mem 

MULTIPLY  BY  2 

RCH 

0A02 

o  INTO  6,  A  I NTS  X 

STB 

MEM 

restore  mem  to  zero 

add 

three 

read  sr  write  switch 

STa 

cntrl 

set  switch 

lda 

TQP12/2 

retrieve 

PRBPeR  tape  OPERATION 

STA 

RWT 

LDa 

RDT,2 

obtain  return  address 

STA 

RBT 

LDX 

ebadr 

ADDRESS  -  ERASABLE  core 

skn 

CNTRL 

STB 

*3*2 

clear  kfile  indicator 

lda 

*  1  *  2 

2ND  ARGUMENT 

sta 

RCDUNT 

WQRD  C0UNT 

etr 

maxwc 

CLEAR  NEGATIVE  WQRD  COUNT 

RSH 

10 

EXTENDED  W0RD  COUNT 

MRG 

TSPil 

EXTENDED  W8RD  C8UNT-E0M 

STA 

MODE 

STB 

PW8RD 

LDa 

*E»ADR 

ADDRESS  -  1ST  ARGUMENT 

etr 

adrmsk 

mrg 

PW0RD 

normal  wqrd  count 

STA 

PWQRD 

INTERPLACE  POT  WORD 

LDA 

*2>  2 

3RD  ARGUMENT 

STA 

magt 

logical  magtape  unit 

MRG 

RWT 

merge  with  tape  operation 

STA 

RWT 

SKN 

CNTRL 

READ  -  WRITE  SWITCH 

BRU 

RERUN 

LDA 

T8P10 

beginning  of  tape  test 

MRG 

magt 

STA 

$  +  2 

BRM 

DELAY 

N8P 

TAPE  OPERATION 

BRM 

ERAS 

BRU 

RERUN 

ERAS 

PZE 

erase  TAPE  forward 

lda 

T8P9 

MRG 

MAST 
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